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Abstract 
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GENERATION OF A DESIRED THREE-DIMENSIONAL ELECTROMAGNETIC 
FIELD 



FIELD OF THE INVENTION 

The present invention relates to a method and a system for synthesizing a prescribed 
three-dimensional electromagnetic field based on generalized phase contrast 
imaging. 

BACKGROUND OF THE INVENTION 

It is well known to form an image by phase contrast imaging methods in which phase 
modulation of light is converted into intensity modulation. As opposed to intensity 
modulation, phase modulation does not involve loss of energy. 

A generalized phase contrast imaging method and system for synthesizing a 
prescribed intensity pattern Is disclosed in WO 96/34207, which is hereby 
incorporated by reference. The generalized method is not based on the so-called 
Zerntke approximation that the phase shift <p is less than 1 radian. An improved 
method is provided without this approximation and based on imaging with a simple 
one-to-one mapping of resolution elements or pixels of a spatial phase modulator 
and resolution elements of the generated Intensity pattern. 

The disclosed phase contrast imaging method of synthesizing a prescribed intensity 
pattern /(*',/), comprises the steps of: 

• dividing the intensify pattern I(.x',y') « lo(*',jOf into pixels in accordance with 
the disposition of resolution elements (x,y) of a spatial phase mask having 

a plurality of individual resolution elements (x,y), each resolution element (x,y) 
modulating the phase of electromagnetic radiation incident upon it with a 
predetermined phasor value e'*^', 

• radiating electromagnetic radiation towards the spatial phase mask, 

• Fourier or Fresnel transforming the modulated electromagnetic radiation, 

• phase shifting with a spatial phase filter (SPF) in a region of spatial frequencies 
comprising DC in the Fourier or Fresnel plane, the modulated electromagnetic 
radiation by a predetermined phase shift value 8 in relation to the remaining part 
of the electromagnetic radiation, and 

• forming the intensity pattern by Fourier or Fresnel transforming, respectively, the 
phase shifted Fourier or Fresnel transformed modulated electromagnetic 



radiation, whereby each resolution element (x,y) of the phase mask is imaged on 
a corresponding resolution element fry) of the image, 

• calculating the phasor values e'* ( *' y) of the phase mask and the phase shift value 
9 in accordance with 

for selected phase shift values 8, « being the average of the phasors e'**> of 
the resolution elements of the phase mask, 

- selecting, for each resolution element, one of two phasor values which represent 
a particular grey level, and 

• supplying the selected phasor values d** to the resolution elements (x,y) of the 
spatial phase mask. 

In one embodiment disclosed in WO 96/34207, the spatial phase mask is positioned 
at the front focal plane of a lens while the spatial phase filter Is positioned in the back 
focal plane of the lens, whereby a first electromagnetic field at the spatial phase 
mask is Fourier transformed by the lens Into a second electromagnetic field at the 
spatial phase filter. Thus, specific spatial frequencies of the first electromagnetic field 
will be transmitted through the spatial phase filter at specific positions of the phase 
filter. For example, the energy of the electromagnetic radiation at zero frequency 
(DC) is modified and transformed onto the intersecting point of the Fourier plane and 
the optical axis of the lens also denoted the zero-order diffraction region by the phase 
filter. 

SUMMARY OF THE INVENTION 

The above-mentioned method operates on a plane Incoming electromagnetic field 
with the aim to generate a two-dimensional intensity pattern. It Is an object of the 
present invention to provide a method and a system synthesizing a prescribed three- 
dimensional electromagnetic field for further flexibility, for example to be able to focus 
light energy in a three-dimensional volume. Such a method and apparatus may be 
utilized in advanced optical micro- and nano-manlpulation, such as by provision of a 
multiple-beam optical tweezer. 

According to a first aspect of the Invention the above-mentioned and other objects 
are fulfilled by a phase contrast system for synthesizing an output electromagnetic 
field u(x", y", z"), comprising 
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a first phase modifying element for phase modulation of an input electromagnetic 
field by phasor values e my \ 

first Fourier or Fresnel optics for Fourier or Fresnel transforming the phase 
modulated electromagnetic field positioned in the propagation path of the phase 
5 modulated field, 

a spatial filter for filtering the Fourier or Fresnel transformed electromagnetic 
radiation by 

in a region of spatial frequencies comprising DC in the Fourier or Fresne! 
plane 

1 0 phase shifting with a predetermined phase shift value 9 the modulated 

electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant 6 t and 

15 in a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the ampIJtude of the modulated electromagnetic radiation 
with a constants, 

second Fourier or Fresnef optics for forming an electromagnetic field o(x', y') by 
Fourier or Fresnel transforming the phase shifted Fourier or Frresnel transformed 
20 electromagnetic field, and 

a second phase modifying element for phase modulating the electromagnetic field 
a(x\ y'} into the output electromagnetic field o(x\ yjeF (iC y,) propagating as the 
desired output electromagnetic field u (x \y'\z") . 

According to a second aspect of the present Invention, the above and other objects 
25 are fulfilled by a method of synthesizing an output electromagnetic field u(x" t y*\ z") t 
comprising the steps of 

phase modulating an input electromagnetic field by phasor values e mv \ 

Fourier or Fresnel transforming the phase modulated electromagnetic field, 

filtering the Fourier or Fresnel transformed electromagnetic radiation by 

30 in a region of spatial frequencies comprising DC in the Fourier or Fresnel 

plane 
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phase shifting with a predetermined phase shift value e the modulated 
electromagnetic radiation In relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
5 with a constant B, and 

in a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant A, 

forming an electromagnetic field o{x', /) by Fourier or Fresnel transforming the phase 
1 0 shifted Fourier or Fresnel transformed electromagnetic field, and 

phase modulating the electromagnetic field o(x', y') Into the output electromagnetic 
field o(x\ y')e w t* '•« propagating as the desired output electromagnetic field 

The method may further comprise the steps of 

15 dividing the electromagnetic field o(x\y') into pixels In accordance with the disposition 
of resolution elements (x, y) of a first phase modifying element having 

a plurality of Individual resolution elements (x, y), each resolution element 
(x, y) modulating the phase of electromagnetic radiation incident upon it with 
a predetermined phasor value B iw *' v, l 

20 calculating the phasor values e lw *- v > of the phase modifying element and the phase 
shift value $ substantially in accordance with 

wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
25 order diffraction region, 

B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region, 



a = 



a 



' ex Pfe) ls the average of the phasors e® m af the resolution elements of the 
phase modifying element, and 

30 0~<£-<k, and 
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10 



K = 1 - J 0 (1 ._22 jt?) • wherein 
Jo Is the zero-order Bessel function, and 

rj relates the radius R, of the zero-order filtering region to the radius R 2 of the main- 
lobe of the Airy function of the input aperture. ?/ - if, IR, m (0.6 l)" 1 ArA/ r , 

selecting, for each resolution element, one of two phasor values which represent a 
particular grey level, and 

supplying the selected phasor values e H «>» to the respective resolution elements 
{x, y) of the first phase modifying element, and 

supplying selected phasor values e w « to respective resolution elements (x\ y*) of a 
second phase modifying element having a plurality of individual resolution elements 
(x\ y'), each resolution element <x\ y') modulating the phase of electromagnetic 
radiation Incident upon it with the respective phasor value for generation of the 
output field Q(x\y')e W[x '-*\ 

The mathematical expressions will be further explained below. 

1 5 The axis of propagation of a plane electromagnetic field is perpendicular to the 
electric and magnetic fields. 

It should be noted that, in each resolution element of the first phase modifying 
element, one of two phasor values which represent a particular grey level of the 
amplitude component of the eiectromagnetic field o(x',y') may be selected. 
20 In an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields incident upon it outside the phase shifting 
regions, I e. A Is equal to one Dr approximately equal to one, 

in an embodiment of the present invention, the spatial phase filter substantially does 
not attenuate the electromagnetic fields Incident upon it inside the phase shifting 
25 region. i.e. Sis equal to one or approximately equal to one. 

It is also preferred that the phase shift value 0 substantially fulfils the equation 



*l*fl 



for a lossless filter with A = 1 and B = l . In a preferred embodiment of the present 
invention, the phase shift 8 is equal to % or approximately equal to * Accordingly the 
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prevlous equation leads to Ka~/ Z and the phase values, M*>y), of the first phase 
modifying element may be calculated in accordance with 

iCA"'JJsin( ( *(x^))&^ = 0 

1 A 

where a Is the illuminated area of the first phase modifying element 

The electromagnetic field or radiation may be of any frequency range of the 
electromagnetic spectrum, i.e. the gamma frequency range, the ultraviolet range, the 
visible range, the infrared range, the far Infrared range, the X-ray range, ihe 
microwave range, the HF (high frequency) range, etc. The present invention is also 
applicable to particle radiation, such as electron radiation, neutron radiation, eta 

Preferably, the electromagnetic fields are monochromatic or quasi-monochromatic so 
that the energy of the electromagnetic fields is concentrated in a narrow frequency 
bandwidth. Since the phase contrast generated amplitude pattern is reconstructed by 
interference of two electromagnetic fields generated by different phase shifting of 
different parts of the incoming field, it is required that the frequency range of the 
emitted electromagnetic field is sufficiently narrow to ensure that the two 
electromagnetic fields are coherent so that their superposition generates the desired 
amplitude pattern. If the frequency range is too broad, the two fields will be 
incoherent and the phase Information will be lost as superposition of non-coherent 
fields results in a summation of fhe intensities of the two fields, it is required that the 
difference between individual delays of electromagnetic fields to be superpositioned 
is less than the wavelength of the fields. This is a relaxed requirement that allows the 
electromagnetic fields to be relatively broad-banded. For example in fhe visible 
range, a Xe-lamp or a Hg-lamp can be used as a light source in a system according 
to the present invention with the advantage compared to a laser light source that 
speckle noise is reduced. The requirements of the spatial coherence of the 
electromagnetic fields depend upon the space bandwidth product of the 
corresponding system and how close the required system performance is to the 
theoretically obtainable performance of the system. 

Preferably, the electromagnetic radiation is generated by a coherent source of 
electromagnetic radiation, such as a laser, a semi-conductor laser, a strained multi- 
quantum well laser, a vertical cavity surface emitting laser (VCSEL), a maser, a 
phase-locked laser diode array, a light emitting diode, a pulsed laser, such as a sub- 
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picosecond laser, etc, or an array of such sources, However, as already mentioned, 
a high-pressure arc lamp, such as an Hg lamp, a Xe lamp, etc, may also be used and 
even an incandescent lamp may be used as a source of electromagnetic radiation. 

Each phase modifying element changes the phase of an electromagnetic field 
5 incident upon it. Optionally, It may also change the amplitude of an electromagnetic 
field incident upon It. Each phase modifying element may transmit or reflect the 
incident electromagnetic field. Each phase modifying element may be divided into a 
number of resolution elements, each of which modulates the Incident electromagnetic 
field by changing its phase by a specific predetermined value. The predetermined 

1 0 values are assigned to each resolution element In different ways depending upon the 
technology applied in the component. For example In spatial light modulators, each 
resolution element may be addressed either optically or electrically. The electrical 
addressing technique resembles the addressing technique of solid-state memories in 
that each resolution element can be addressed through electronic circuitry to receive 

1 5 a control signal corresponding to the phase change to be generated by the 

addressed resolution element. The optical addressing technique addresses each 
resolution element by pointing a light beam on it, the intensity of the light beam 
corresponding to the phase change to be generated by the resolution element 
illuminated by the light beam. ' 

20 Spatial phase modulation may be realized utilizing a fixed phase mask, a liquid 
crystal device based on liquid crystal display technology, a MEMS (micro electro- 
mechanical system), a MOEMS (micro opta-electro-mechanical system), such as a 
dynamic mirror device, a digital micro-mirror array, a deformable mirror device, etc, a 
membrane spatial light modulator, a laser diode array (integrated light source and 

25 phase modulator), smart pixel arrays, etc. 

Seiko-Epson produces a transmitting liquid crystal SLM (LC-SLM) having a high 
resolution matrix of transparent liquid crystal elements wherein the relative 
permittivity of each element can be electrically modulated in order to vary the 
refractive index and thereby the optical path length of the element. 

30 Meadowlark produces a parallel-aligned liquid crystal (PAL-SLM) with a high fill 
factor, but this device has a very low resolution in that it contains only 137 phase- 
modulating elements. 

Hamamatsu Photonics produces a dynamically controllable PAL-SLM with VGA or 
XGA resolution. 
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Texas Instruments produces a Digital Mirror Device (DMD) having an array of mirrors 
each of which can be tilted between two positions. 

The spatial phase filter is typically a fixed phase mask, such as an optically flat glass 
plate coated with a dielectric layer in the region wherein the modulated 
5 electromagnetic field is phase shifted 9 in relation to the remaining part of the 
electromagnetic field. However, the spatial phase modulators mentioned in the 
previous section may also be used for spatial phase filters. In addition, non-linear 
materials providing self-phase modulation, such as Kerr-type materials, can also be 
used for introducing the phase shift 9. 

1 0 An imaging system maps the phase modulating resolution elements of the first phase 
modifying element onto the second phase modifying element. This imaging system 
may comprise a 4f-lens configuration (two Fourier transforming lenses utilizing 
transmission of light or one Fourier transforming lens utilizing reflection of light) or a 
single imaging lens. However, any optical imaging system providing a filtering plane 

1 5 for the spatial phase filter may be applied in a phase contrast imaging system. 

In the method and system according to the present Invention, the electromagnetic 
field o(x', y) is generated by superposition of electromagnetic fields in the Image 
plane of the imaging system. The first phase modifying element changes the phase 
values of an electromagnetic field incident upon it and the imaging system directs the 

20 electromagnetic field with changed phases reflected from or transmitted through the 
phase modifying element towards the spatial phase fliter. The phase filter phase 
shifts a part of the electromagnetic field and the imaging system is adapted to 
superimpose in the image plana the phase shifted part of the electromagnetic field 
with the part of the electromagnetic field that is not phase shifted by the spatial phase 

25 filter. 

According to a preferred embodiment of the invention, the first phase modifying 
element is positioned at the front focal plane of a lens while the spatial phase filter is 
positioned in the back focal plane of the lens, whereby a first electromagnetic field at 
the phase modifying element is Fourier transformed by the lens Into a second 

30 electromagnetic field at the phase filter. Thus, specific spatial frequencies of the first 
electromagnetic field will be transmitted through the spatial phase filter at specific 
positions of the phase filter. For instance, the energy of the electromagnetic field at 
zero frequency (DC) is transmitted through the phase filter at the intersecting point of 
the Fourier plane and the optical axis of the lens also denoted the zero-order 

35 diffraction region. 
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It is an advantage of the invention that utilisation of arrays of sources is facilitated in 
that the positioning and/or shapes of the phase shifting regions of the phase filter 
may be matched to the geometry of the source. For example, if a linear array of 
VCSELs forms the source, the phase shifting regions of the spatial phase filter form a 
5 corresponding linear array of phase shifting regions, each of the regions being 
positioned at the zero-order diffraction region of a respective VCSEL in the VCSEL 
array. Further, the shape of each phase shifting region may match the shape of the 
zero-order diffraction region of the respective VCSEL 

Likewise, a phase filter may match a source with a specific geometrical shape with a 
1 0 continuous phase shifting region covering an area of the phase filter that corresponds 
to the zero-order diffraction region of the source. 

Thus, the energy of the electromagnetic fields of the system may be distributed over 
a large area compared to the area of a zero-order diffraction region of a single plane 
electromagnetic field of a known phase contrast Imaging system. 

1 5 Thus, the phase shifting regions of the spatial phase filter may form a rectangular 
array, a circular array, a linear array, two linear crossing arrays, a continuous region, 
a ring, etc. 

At least two substantially plane electromagnetic fields with different axes of 
propagation may be generated in a time multiplexed manner, e.g. by a scanning 
20 mirror or prism, deflecting or reflecting a beam of electromagnetic field In different 
directions of propagation, 

The capability of handling high energy levels of electromagnetic fields of the present 
Invention may be utilized for provision of a 3D laser cutter comprising a system 
according to the present invention. 

25 Further, the capability of handling high energy levels In combination with the 
capability of generating a desired three-dimensional field comprising desired light 
beams may be utilized for provision of an optical tweezer or an array of optical 
tweezers according to the present Invention. 

In an embodiment of the present invention, wherein the apertures of the system is of 
30 insignificant Importance to the operation of the system and calculation of the phasor 
values, K is equal to one or approximately equal to one. 

For a more detailed understanding of the invention, the Zemike approximation is 
reviewed below, followed by a generalization where the above-mentioned 



-10- 



mathematica! expressions are derived for an on-axis centred phase contrast filtering 
implementation. 

The Zernike phase contrast method aliows for the visualization of phase 
perturbations by the use of a Fourier plane phase shifting filter. The Dutch physicist 
5 Fritz Zernike received the Nobel Prize in 1 953 for inventing this method, which led to 
a break-through in medicine and biology by making essentially transparent cell or 
bacteria samples clearly visible under a microscope. Its successful operation, 
however, requires that the spatial phase distribution. $(x,y) , at the input is limited to 
a "small-scale" phase approximation where the largest phase is typically taken to be 
1 0 significantly less than nl 3 .According to this assumption, a Taylor expansion to first 
order is sufficient for the mathematical treatment so that the Input wavefront can be 
written as 

The light corresponding to the two terms in this "smalt-scale" phase approximation 
1 5 can be separated spatially by use of a single lens where the phase distribution is 
located in the front focal plane and the corresponding spatial Fourier transformation 
is generated in the back focal plane of the lens. In this first order approximation the 
constant term represents the amplitude of on-axis light focused by the lens In the 
back focal plane and the second spatially varying term represents the off-axis light. 
20 Zernike realised that a small phaBe shifting quarter wave plate acting on the focused 
light makes It possible to obtain an approximately linear visualization of small phase 
structures by generating interference between the two-phase quadrature terms In Eq. 
(1): 

/(*',j' , )*I + 2 l |>(* , ,y) {2) 

25 It should be noted that a three-quarter waveplate works equally well to produce 
contrast, but the plus sign in Eq. (2) is negated leading to so-called negative phase 
contrast. A substantial improvement in the visibility of the Zernike phase contrast 
visualization in Eq. (2) requires strong damping of the focused light In addition to the 
phase shift required to generate the contrast. 

30 In the general case, where we are not limited to a small-scale input phase 
perturbation we cannot assume that a series expansion to first order as in the 
Zernike approximation Is a sufficient representation of a given phase perturbation. 
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Hi0her order terms in the expansion need to be taken into account, so the expansion 
takes the form: 

However, here the spatially varying terms can not be considered as separate from 
the supposedly focused light represented by the first term in this Taylor series 
expansion, as Is Implied by the Zernike approach, and all of these spatially varying 
terms contribute to the intensity of the on-axis focused light. For a significant 
modulation In the input phase, this contribution of the spatially varying terms can 
result In a significant modulation of the focal spot amplitude in the back focal plane of 
the lens. These terms can in fact result in either constructive or destructive 
Interference with the on-axis light, although the net result will be an attenuation of the 
focused light amplitude, which only has a maximum value for a perfect unperturbed 
plane wave at the input. 

For phase objects not fulfilling the Zernike approximation we must, therefore, find an 
alternative mathematical approach to that of the Taylor expansion given in Eq. (3). 
We have chosen a Fourier analysis as a more suitable technique for completely 
separating the on-axis and higher spatial frequency components. This gives the 
following form for cxp^jc,^)) , where (x,y) eQ: 

exp(/<j>(,*,y)) = JJ^ |fexp(z(Ji(x,7))^+"higher frequency terras" (4) 
U ; a 

In this Fourier decomposition the first term Is a complex valued constant linked to the 
on-axis focused light from a phase object defined within the spatial region, £2, and 
the second term describes light scattered by spatially varying structures in the phase 
object. Comparing Eq. (3) and Eq, (4) It is apparent that the first term of Eq, (3) is a 
poor approximation to the first term of Eq. (4) when operating beyond the Zernike 
small-scale phase reglma. 

An important issue to consider when analysing the effect of spatial filtering of the light 
diffracted by phase perturbations is the definition of what spatially constitutes focused 
and scattered light. In the previous description of Zernike phase contrast it was 
assumed that the focused light is spatially confined to a somewhat unphysical delta 
function. As we know, any aperture truncation inherent In any practical optical system 
will lead to a corresponding spatial broadening of the focused light. It Is therefore 
essential that we define the terms "focused light" and "scattered light" explicitly for 
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such a system. In this context it is necessary to look more carefully at the sequence 
of apertures confining the lightwave propagation through a typical optica! set-up. 

A commonly applied architecture that provides an efficient platform for spatial filtering 
is illustrated In Fig. 1 and is based on the so-called 44 configuration.. An output 
interferogram of an unknown phase object or phase disturbance is obtained by 
applying a truncated on-axis filtering operation in the spatial frequency domain 
between two Fourier transforming lenses (L1 and L2). The first lens performs a 
spatial Fourier transform so that directly propagated light Is focused into the on»axis 
filtering region whereas spatially varying object information generates light scattered 
to locations outside this central region. We can describe a general Fourier filter in 
which different phase shifts and amplitude damping factors are applied to the 
"focused" and "scattered* light In Fig 1, we show a circularly symmetric Fourier filter 
described by the amplitude transmission factors A and B for the "scattered" and 
"focused" light respectively and by the relative phase shift ft These filter parameters 
can be chosen to replicate any one of a large number of commonly used filter types 
(he. phase contrast dark central ground, point diffraction and field-absorption 
filtering). By applying a given Fourier filter and a second Fourier lens, we obtain an 
interference pattern In the observation plane. The focused on-axis light acts as the 
synthetic reference wave (SRW) In the common path Interferometer (CPI) system, 
this interferes with the scattered light to generate the output Interference pattern, In 
the following section we discuss the Importance of the SRW and show how it 
influences, among other things, the choice of the Fourier filter parameters. 

Having described the generic optical system that makes up the CPJ, we turn to a 
detailed analytical treatment of the important elements in this system. Assuming a 
circular input aperture with radius, Ar , truncating the phase disturbance modulated 
onto a collimated, unit amplitude, monochromatic field of wavelength A.we can 
describe the incoming light amplitude a(x>y) by, 

a (xy)=oirc(r/&r ) exp (fy (x, y )) (5) 

at the entrance plane of the optical system shown In Fig, 1 using the definition that 
the circ-function is unity within the region, r -^x l + y l <Ar , and zero elsewhere. 
Similarly, we assume a circular on-axis centred spatial filter of the form: 



(6) 
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where B <=[(); l] is the chosen filter transmittance of the focused light, 9 e[0;2/r] is 
the applied phase shift to the focused light and A e[0;l] Is a filter parameter 
describing field transmittance for off-axis scattered light as indicated in Fig. 1 The 
spatial frequency coordinates are related to spatial coordinates in the filter plane 

5 such that; (/, ./„) - WV (*, ,y,) and f, = ffiTjj . 

Performing an optical Fourier transform of the input field from Eq. (5) followed by a 
multiplication with the filter parameters in Eq. (6) and a second optical Fourier 
transform (corresponding to an inverse Fourier transform with inverted coordinates) 

we obtain an expression for the intensity /fV J y , )=|c?(*', < y , )| 3 describing the 

1 0 interferogram at the observation plane of the 44 set-up: 

where g(r') is the synthetic reference wave (SRW) and the terms a and $(x',y') 
are given by: 
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It should be noted that to achieve a tractable analytic expression In Eq. (7) it has 
been assumed that the spatial frequency content of the phase object is sufficiently 
described by the term, a , within the on-axis centred filtering region characterized by 
20 the spatial frequency range Af r . 

The generally complex valued and object dependent term, a , corresponding to the 
amplitude of the focused light plays a significant role in the expression for the 
interference pattern described by Eq. (7). Referring to the discussion in the 
introduction, we are now able to confirm that the frequent assumption, that the 
25 amplitude of the focused light Is approximately equal to the first term of the Taylor 
expansion in Eq. (1), can generally result in misleading Interpretations of the 
interferograms generated at the CPI output. 

Of similar Importance in the analysis of Eq. (7) is the term g(r') describing the spatial 
profile of the SRW, diffracted from the aperture formed by the on-axis centred filtering 
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region. It is the interference between this SRW term, carrying the information about 
the filtering parameters, and the Imaged phase object that generates the output 
interferogram. Thus, it is important to obtain an accurate description for the SRW and 
thereby an accurate derivation for Eq. (7). The zero-order Hankei transform followed 
by a series expansion in the spatial dimension, f, wifl be used to describe the SRW. 
This is a relatively simple approach, which to the knowledge of the author has not 
previously been applied to this problem. 

For a circular Input aperture with radius, Ar , we can describe the radius of the 
corresponding central phase shifting region of the Fourier filter (characterized by the 
parameters 6 and in terms of a radial spatial frequency range &f r . We can thus 

obtain the following expression for tie SRW by use of the zero-order Hankei 
transform: 

g(r'h^fM 2 * A rty«(2nr<f r )df, (9) 

In order to simplify the analysis, we introduce a term which explicitly relates the 
radius of the central filtering region, 2?, , to the radius of the main-lobe of the Airy 
function, R, , resulting from the Fourier transform of the circular input aperture alone. 
We can thus express 7 In terms of At and Af f such that: 

T\ = RjR l ~(0.6if , &rAf r (10) 

where the factor of 0.61 arises from the radial distance to the first zero crossing of 
the Airy function corresponding to half of the Airy mainlobe factor, of 122. If we make 
this substitution in Eq.(0) and then perform a series expansion In r' , we obtain the 
following expression for the SRW: 

g{r ') «= W 0 (l.22«ti)-[(0.61nn) S h (l&m\)](r '/Arf + 
j[(0.6l7tT ) ) 3 /4][2J 3 (lJ2K n )~0-61nnJ- 4 (l.22KT } }]J(r'/A^) 4 .„, 

In this expansion, the SRW is expressed in radial coordinates normalised to the 
radius of the imaged input aperture. This can easily be scaled to allow for a 
magnification within the imaging system, though for the remainder of our analysis a 
direct imaging operation is assumed. From Eq.(11 ) it Is apparent that the SRW will 
change as a function of the radius of the central filtering region. Additionally, it is 
clear that the SRW profile is not necessarily flat over the system output aperture. 



This is an important, yet often neglected, factor in determining the performance of a 



Depending on the accuracy needed for the description of the interferograms one can 
choose to Include a number of spatial higher order terms from the expansion in Eq, 
5 (11). The influence of the higher order terms has the largest impact along the 
boundaries of the imaged aperture. For q -values smaller than 0.627 and when 
operating within the central region of the image plane, spatial higher order terms are 
of much less significance and we can safely approximate the synthetic reference 
wave with the first and space invariant term in Eq. (1 1 ): 

10 g(r'ectmtral region) xl~J 0 (l22nx}) (12) 

so that we can simplify Eq. (7) to give: 



where K= l~./ 0 (L22«rj;), The influence of the finite on-axis filtering radius on the 
focused light Is thus effectively included as an extra "filtering parameter" so that the 
1 5 four-parameter filter set (A, B, 6, K(ij)) together with the complex object dependent 

term, a , effectively defines the type of filtering scheme we are applying. 

Having determined a suitable operating range for the CPI in terms of the production 
of a good SRW, we must now examine the role that the remaining filter parameters 
play in the optimisation of a CPi. 

20 From Eq. (1 3) we see that the filter parameters (A, B, $) can be combined to form a 
single complex valued term, C, the combined filter term, such that: 



CPI. 



l{x\f) s A^expffifayty+K^BA-* cxp(i8)-l)[ 2 



(13) 



C~\C\ exp (ty c ) = BA' X exp (&) - 1 



(14) 



therefore, Eq. (13) can be simplified to give: 



(15) 



25 where 



BA' 1 = ^/l+2|C|cos(^ c )+|C| 2 
sk0 = (^-')"'|C|sm(^ c ) 



(16) 
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Since it is a complex variable* the combined filter term C r which effectively describes 
the complex filter space, can be considered to consist of a vector of phase y/ c and 

length |Cj as shown in Eq,{14)~ Thus in order to obtain an overview of the operating 
space covered by all the possible combinations of three independent filter 
5 parameters (A>B>0) we can now instead choose to consider a given filter in terms of 

the two combined parameters and [cj. However, referring to Eq,{15), it can be 
seen that the fitter parameter, A, also appears independently of the combined filter 
term, G. Fortunately, this issue can be resolved by considering that the term SA ~ l 
from Eq. (14) must be constrained in the following way: 

10 

BA~ l <l r> A = \ fl=|C + 1] 
. si -± A = ), B«1 (17) 

These constraints arise from the adoption of a maximum irradlance criterion 
minlmisfng unnecessary absorption of light in the Fourier filter, which reduces both 
1 5 irradiance and the signal to noise ratio in the CPI output 

In the previous sections we derived expressions relating the spatial average value of 
a given phase disturbance to obtain peak irradlance and optimal visibility in 
combination with high accuracy in systems with unknown wavefront phase 
disturbances. We saw that if a CPI is applied to wavefront sensing or the 

20 visualisation of unknown phase objects the Generalised Phase Contrast (GPC) 
method specifies the filter phase and aperture size parameters for achieving optimal 
performance in extracting and displaying the phase information carried by the 
incoming wavefront On the other hand, in cases where we have control over the 
incoming wavefront or phase modulation the GPC method provides extra means of 

25 optimisation by encoding the phase distribution itself in addition to modifying the filter 
parameters. The two main scenarios: A) synthesizing the spatial phase for intensity 
display or B) measuring the spatial phase with high accuracy, strongly influences 
which of the parameters in the analysis that should be kept fixed and which could be 
changed or adapted, The first approach is particularly useful when the filter 

30 parameters have a restricted dynamic range or are fixed* The rigorous derivation of 
the equations for choosing these parameters will be derived in this section* 

When synthesizing an input phase distribution for optimal visibility of an output 
intensity pattern the situation is more relaxed than the situation involving accurate 
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interferometrlc measurements of unknown phase disturbances, The parameter r\ can 
therefore in most cases be chosen to completely encompass the zero-order light with 
the result that the term, K , tends to unity as the Bessel function tends to zero in Eq. 
(12), For this particular case, the SRW becomes a flat top profile where we can 
5 achieve nearly 1 00% light efficiency. For smaller and irregular phase patterns fine- 
tuning of ?] in the region 0.4-0,6 provides for an efficient operating regime while 
maintaining minimal losses. 

in order to optimize a synthesized light distribution for maximum contrast, we wish to 
generate an intensity distribution with a lowest Intensity equal to zero, i,e. in at least 
10 one point [x^y f Q ): 

/(^.34i*o) = 0 (20) 

where fa is the relative phase shift generating zero-intensity in (aj.^oftfie 
obsen/ation plane. 

Applying this dark background condition In Eq, (13) we can obtain the following 
1 5 expression for a no-loss phase-only filter with filter transmission 
parameters,/! = 5*4: 

^p|(l-exp(i8))-cxp(4) (21) 

A key point arising from Eq. (21) is that we now have a simple way of expressing a 
new design criterion relating the spatial average value of any input phase pattern to 
20 the zero-order phase shift of a matched Fourier phase filter. Since K is by definition 
positive and by taking the modulus of Eq. (21 ) we obtain; 

i^|a|=|2sm(e/2)P (22) 

Eq. (22) is a key result for the fully transmisslve wavefront engineered GPC mapping 
that makes it possible to deduce the range of valid phase parameters fulfilling our 

25 design criteria from Eq. (20). The largest possible value that the term, £|<£| , takes on 

is unity, this leads to the following solution interval for Eq. (22) within a full phase- 
cycle; 

8 = [7t/3;57t/3] (23) 
From Eq. (22) we also observe that id#| can take on a value limited to the interval: 
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jr|5|-lK { i]'. < 24 > 

Eq. (22) and the solution intervals described by Eqs. (23)-(24) specify the design 
parameters for achieving optimal performance in extracting and displaying the phase 
information carried by the Incoming wavefront. Moreover, Eq, (22) hints towards extra 
5 means of optimisation by encoding the phase modulation depth itself in addition to 
the no-loss phase-only filter. This last approach is particularly useful when the filter 
phase has a restricted dynamic range or is fixed. 

Now, assuming that we have a fixed and fully transmissive phase-only filter, the best 
choice for the filter parameter is a value that allows for the largest dynamic range of 
1 0 phasor values at the input Accordingly, the smallest possible real value, Ka = I / 2 , 
is desirable implying that 9 = k , leading to the output intensity distribution: 

/(x\/) = 2[t-cos(H*',/))] (25) 

Inserting Ka = 1 / 2 in Eq. (8) we obtain the following two requirements for the Input 
encoded phase function <j>{x,y) : 
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KK A Jjcos(0f>,y))fo^ m X 

(26) 



where A is the Illuminated area of the phase modifying element. 

We observe that It is only the first requirement In Eq. (26) that is directly related to the 
output intensity in Eq. (25) via the cosine term. Since there are always two choices 
for a given phasor value that result in the same cosine value (excluding 0 and n), we 
20 notice that the second requirement can subsequently be fulfilled independently of the 
first requirement simply by complex conjugating an appropriate number of phasor 
values. This feet is a key feature of the GPC-method since it makes it possible to 
solely concentrate on the first requirement In the process of synthesizing a desired 
an virtually no-loss grey level intensity pattern. 

25 The first requirement in Eq. (26) can be fulfilled by several means, Including: dynamic 
phase range adjustment, fill factor encoding, phase-histogram adjustment, spatial 
scaling of phasor pattern, raster encoding etc. in a histogram adjustment technique 

one will typically start out with a desired relative Intensity distribution /(a'j/L,,^ 

where the maximum achievable intensity level is unknown but relative intensity levels 
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are known and the lowest intensity level is fixed by the background criterion of Eq. 

(20), The procedure is now to adjust the histogram for /{*',/) white 

maintaining identical relative Intensity level ratios until the first requirement in Eq, (26) 
Is fulfilled. Subsequently, the second requirement In Eq. (26) is fulfilled by complex 
conjugating an appropriate part of the phasors. The simplest procedure is to complex 
conjugate every second identical phasor value independently of the spatial location, 
However, this "phasor flipping 11 procedure can also be turned into an advantageous 
tool (an extra degree of freedom) for manipulating the spatial frequency content in 
order to optimize the separation of low and high spatial frequency terms at the 
Fourier filter plane by taking the spatial phasor location into account Eg, 
neighbouring phasor values can be chosen to have a maximum difference between 
them, thereby introducing high spatial frequency modulation easing the filtering in the 
spatial Fourier domain. 

In most cases, however, equalized output intensity levels are sufficient In the 
succeeding analysis, we therefore focus on the encoding of the Input phase levels to 
achieve binary output intensity levels, A derivation based on ternary phase levels 
allows for the widest range of binary intensity pattern encoding and automatically 
provides for the simplified but important binary phase level encoding as a special 
case. 

For the ternary phase encoding, we consider the illuminated portion of the input 
aperture area, A f as divided into sub-areas A 0i Ai and A 2 with respective phase 
values $ 0 , A and 0 2 , We are aiming for the derivation of general expressions 
relating the addressing parameters for the phase modulation to the range of possible 
phase parameters of the Fourier filter obeying the design criterion we have already 
set out We can express the total truncated area and its average phase modulation, 
as the sum of the phase-weighted sub-areas: 

A 0 exp(4) + A, oq>(tf ) + A 3 exp(4) * a|o| (27) 

This can be further simplified by expressing the sub-areas as fractions of the total 
area, A such that F } -A } /A and = A 2 / A ; 

{\-P x -F 2 ) exp (4 ) + F { exp ) + F 2 exp (% ) « |a| (28) 

As previously mentioned we are interested in binary intensity patterns with levels 
corresponding to the input phase values. In this case the dark background region is 
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10 



defined by (a 0 ,#>) and the bright output level of intensity, /« Is determined by 

(A,,$) and (A. 2 ,$,) at the input plane. For the binary output intensity condition it 
follows that: 



(29) 



This equality corresponds to a symmetric condition that can be easily verified 
applying the phasor chart analysis technique demonstrated in section 4. Due to this 
symmetry we can simplify the analysis by applying the following substitution: 

A<i> = $,-$„ =| 0 -Ia (30) 

so that Eq.{28) can be rewritten as: 

F, {swimytf+Fifapi-i^)-!)**-* (l-exp(/9))"' -1 (31) 

it is now a straightforward task to solve Eq. (31) for the real part and the Imaginary 
part respectively, to obtain the following set of equations: 
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F, + F 2 = (2K-1)(2K(]-qos(A$ 
F, -F 2 =sin(0)(2,S:si!i(A(|))(l-cos(e)))" 1 



This can also be expressed in terms of the fractional areas, such that: 



F, ! =(4A'■)"'[(2^~l)(l--coB(A^.))" , +sin(0)(sra(A^)(l-cos(8)))" 1 ] 
F l =(4^)" , [(2i?-l)(i-cos(A^))' t -dn(0)(sin(A<^))(l-cos(^)))" , ] 



(32) 
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20 Since we have focused on solutions where Identical Intensity levels are obtained in 
both the F r region and the F r region We can define the resulting illumination 
compression factor, C, in the following way: 



C - (F, + F 2 = (l - (2Kf )" (1 - cos (A*)) 



(34) 



The minimum compression factor corresponds to uniform illumination at the output 
25 such that Fj + F a = 1 , whereas the maximum compression factor is found to be 
C-s-wfor/sr*^. 
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An interesting special case can be deduced from Eq„ (31) by setting F 2 « 0 , where 
we find that: 

F = F X « (it (l -exp(£8» -l)(JS: (l - wp - exp (»))) ' (35) 
implying that for the binary phase modulation case we must have: 
5 A<H0 (36) 

in order for the flfl factor, F, to be reai-valued. 

This result turns out to be the special case that corresponds to the set of solutions 
where a binary phase pattern serves as the input, 

The second phase modifying element imposes a spatial phase shift, y{x\y'}, on the 
1 0 phase contrast generated electromagnetic field o(x\ y'). This superposition generates 
an arbitrary controllable complex field, o{x\y*)e W)f) with amplitude given by the 

square-rooted Intensity, \o(x\y)\-^jl(x\y t ) i and phasor given by the remaining 

part under division with the amplitude term, o(x\y)e ¥l ^ y<) /\o(x l ,y r )\ - The resulting 

arbitrary controllable complex field, o {x\ y , )e^ x ' t} ^ is capable of rs-dlstributing the 

15 light into an arbitrary three-dimensional focusing within a selected volume of 

operation, This focusing within a volume Is a result of complex wave propagation and 
can be deduced by use of Maxwell's equations. A simplified scalar description of this 
complex wave propagation into a 3D field distribution can be obtained by use of a 
simple piane wave Fourier decomposition: 

20 JJ * ' r (37) 

operating on the Fourier transform, 3, of the controllable complex electromagnetic 
field leaving the second phase modifying element Any subsequent focusing optics 
can be Included in the Fourier decomposition of Eq. (37), 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates schematically a known 4f phase contrast imaging system, 

Fig. 2 illustrates schematically an embodiment of the present Invention , and 

Fig. 3 Illustrates schematically an embodiment of the present invention generating 
5 the desired three-dimentsional field by focusing of o(x f , y')d v (x '- V 'K 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 shows a known 4f CPI phase contrast imaging system 1 . A laser (not shown) 
emits a light beam, which is expanded by a beam expander (not shown) into a plane 
light wave of uniform Intensify and directs It towards a first phase modifying element 

1 0 4. The light beam is transmitted through the first phase modifying element 4 and a 
Fourier transforming lens L1. The first phase modifying element 4 is positioned in the 
front focal plane of the lens L1 and a spatial phase filter 6 }s positioned in the back 
focal plane of the lens L1 that Is also the front focal plane of a lens L2. The Fourier 
transforming lenses L1 , L2 need not have identical focal lengths. Different focal 

1 5 lengths lead to a magnification ratio different from one. The spatial phase filter 6 
phase shifts by 0, and optionally attenuates (by a factor B), the zero order diffraction 
part 8 of the light phase modulated by the first phase modifying element- Optionally, 
the remaining diffraction part of the light modulated by the phase modifying element 
may be attenuated by a factor A. The electromagnetic field o(x', /) is generated in 

20 the back focal plane 9 of the lens L2. 

Figs. 2 and 3 illustrate schematically an embodiment 10 of the present invention, 
comprising a 4f CPI phase contrast imaging configuration as illustrated in Fig. 1 . 
Corresponding parts in the figures are designated with Identical reference numerals. 
It is obvious for the person skilled in the art that the 4f configuration may be 

25 substituted by the 2f or 1f configurations disclosed in WO 96/34207. The operation of 
the 4f CPI phase contrast imaging system is explained with reference to Fig. 1 and is 
not repeated. Again the electromagnetic field o(x\ y') is generated in the back focal 
plane 9 of the lens L2 in front of the second phase modifying element 14. The first 
phase modifying element 4 has a plurality of individual resolution elements (x, y), 

30 each resolution element (x, y) modulating the phase of electromagnetic radiation 
incident upon it with a predetermined phasor value e' M *' y) . As already disclosed, the 
phasor values e^ of the phase modifying element and the phase shift value & are 
calculated substantially in accordance with 



o(^'^')=^[ e xp(^(^'^ , ))- { - i5: R(^" lex P(^)- 1 )] 
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wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction region, 

8 is an optional amplitude modulation of the spatial phase filter in the zero-order 
5 diffraction region, 

a = |«|exp{/^ ) is the average of the phasors e t6( * ,y) of the resolution elements of the 
phase modifying element, and 

$ 88 0 ~ ■ and 

K = 1 - y 0 (1 ,22^7/) , wherein 

10 Jo is the zero-order Bessel function, and 

?1 relates the radius R 1 of the zero-order filtering region to the radius R 3 of the main- 
lobe of the Airy function of the input aperture, ij-RjIt^^ (0.6\y ] &r&f r , 

For each resolution element (x, y), one phasor value is selected from two phasar 
values representing a particular grey level, 

1 5 The computer 1 2 supplies the selected phasor values e l6(x,y) to the respective 

resolution elements (x r y) of the first phase modifying element 4 f and supplies the 
determined 9 value to the spatial phase filter 6, 

The second phase modifying element 14 also has a plurality of individual resolution 
elements (x\ y'), each resolution element (x\ y') modulating the phase of 

20 electromagnetic radiation incident upon it with a predetermined phasor value d WtV) , 
The computer 12 supplies determined phasor values e M *' y l to respective resolution 
elements (x\ y') of the second phase modifying element 14 for modulation of the 
phase of the electromagnetic field o(x' t y 1 ) incident upon it with the respective phasor 
value e M *' ,yP) for generation of the desired output field o(x\ y')e^ < x * rt emitted from the 

25 second phase modifying element 14. 

Thus, any desired amplitude as a function of (x\ y') and any desired phase as a 
function of (x\ /) of the output field o(x\ y')e 141 /) may be synthesized by the 
method and apparatus according to the present invention. 

The output field o(x' p y'Je^^^ propagates and generates the desired three- 
30 dimensional field u(x", y", z"), Le, the light is re-distributed into an arbitrary three- 
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dimensional field distribution within a selected volume. A simplified scalar description 
of this complex wave propagation into a three-dimensional field distribution can be 
obtained by use of a simple plane wave Fourier decomposition: 

where 3 is the Fourier transformation of the controllable complex field o(x\ y'Je" 1 ' !,t '- y,, 
leaving the second phase modifying element 14. 

A further optical system, such as a lens 16, a microscope objective lens, a curved 
mirror, an aspherlc lens, etc, may focus the electromagnetic field afo'.yy {x '" n within 
the selected volume. Any subsequent focusing optics can be Included in the Fourier 
decomposition of the above equation. 

Further, the computer 12 may comprise light control means for controlling the power 
of the light-emitting source generating the field incident on the first phase modifying 
element 4. The computer may also comprise input means, such as a keyboard, a 
mouse, a 3D mouse, 3D virtual reality equipment, a diskette drive, a USB interface, 
an optical disc drive, a network interface, a modem, etc, for receiving a three- 
dimensional field to be synthesized by the system 10. From the received three- 
dimensional field, the computer may be adapted to calculate phasor values e' 4(xy) and 
e ,v *rtfo be transmitted to the respective resolution elements (x, y) and (x\ /) of the 
first and second phase modifying element 4, 14, and the phase shift 9 of the spatial 
phase filter 6 for transmission to the spatial phase filter 6 In accordance with the 
above-mentioned equations. 

PROPOSED APPLICATIONS; 

• 3D optical micro- and nanomanipulation in real-time, 

• 3D optical multi-beam tweezing for manipulation of micro-objects, such as micro- 
components, biological cells, etc, using electromagnetic gradient forces 
proportional to the optical intensity pointing In the direction of the Intensity 
gradient. 

• Optical fractionation, sorting, sifting etc. 

• Efficient and dynamic spot-array generators to provide bias or holding beams for 
3D arrays of photonic elements, such as bistable elements, photonic switches 
and smart pixels. 



• Generation of structured light for machine vision applications. E.g. periodic and 
skew periodic mesh grid illumination in 3D that can be updated in parallel 

• Photolithographic applications (laser 3D direct writing in parallel Without the need 
for sequential scanning), E*g. high power laser direct writing of waveguides in 

5 Ge-doped silica. 

« Volume light intensity modulation in general by use of pure phase modulation 
(radiation focusators). 

• 3D laser beam shaping In real time, 

• 3D image projection without the need for a laser-scanning device, 
10 • Dynamic Infrared Scene Projection (DIRSP). 

• Exposure device for grating and mask production. 

• LIDAR applications. 

• Laser printing in parallel. 
» Laser show applications. 

15 • Atmosphere research. 
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CLAIMS 

1. A phase contrast system for synthesizing an output electromagnetic field 
u(x'\ y", z M ), comprising 

a first phase modifying element for phase modulation of an input electromagnetic 
5 field by phasor values 

first Fourier or Fresnel optics for Fourier or Fresnel transforming the phase 
modulated electromagnetic field positioned in the propagation path of the phase 
modulated field, 

a spatial filter for filtering the Fourier or Fresnel transformed electromagnetic 
10 radiation by 

in a region of spatial frequencies comprising DC in the Fourier or Fresnel 
plane 

phase shifting with a predetermined phase shift value $ the modulated 
electromagnetic radiation in relation to the remaining part of the 
1 5 electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant B r and 

In a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
20 with a constants, 

second Fourier or Fresnel optics for forming an electromagnetic field o(x r , /) by 
Fourier or Fresnel transforming the phase shifted Fourier or Fresnel transformed 
electromagnetic field, and 

a second phase modifying element for phase modulating the electromagnetic field 
25 o(x\ y') Into the electromagnetic field q{x\ y')sF ^ y>) propagating as the desired 
output electromagnetic field u (x n ,y\z") , 

2. A phase contrast system according to claim 1 f wherein 

A*1. 

3. A phase contrast system according to claim 1 or 2, wherein 
30 B = 1. 

4. A phase contrast system according to any of the previous claims, wherein 
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0-17. 

5. A phase contrast system according to any of the previous claims, wherein the 
phasor values e' M * y) of the first phase modifying element and the phase shift value 9 
substantially fulfill that 

o{x\f) s a[sxv($ {x\yty+K$[BA* exp(/0)-l)] 

wherein 

A is an optional amplitude modulation of the spatia! phase filter outside the zero- 
order diffraction region, 

B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region, 

a - |a| exp (z^ ) is the average of the phasors e'*^ y) of the resolution elements of the 
phase modrfying element, and 
^ = $ - $h , and 

K - J - J n (l.22mi) , wherein 

J 0 Is the zero-order Bessel function and 

?? relates the radius ft, of the zero-order filtering region to the radius R z of the main- 
lobe of the Airy function of the input aperture, 7 = R l JR i =* (0-61)" 1 Ar&f r . 

6. A phase contrast system according to any of the previous claims, wherein the 
phase shift value 9 substantially fulfills the equation 

1 1 2|sin0/2| 

7. A phase contrast system according to any of the previous claims, wherein at least 
one of the first and second phase modifying element comprises a complex spatial 
electromagnetic field modulator that is positioned in the path of the input 
electromagnetic field and comprises modulator resolution elements (x m , y m ), each 
modulator resolution element (x mi y m ) modulating the phase and the amplitude of the 
electromagnetic field incident upon it with a predetermined complex value 

am(x m , y m )ew«. 
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8. A phase contrast system according to any of the previous claims, further 
comprising a light source for emission of the Input electromagnetic field, the light 
source comprising a laser array, such as a VCSEL array, 

9. An optical micro-manipulation or multi-beam optical tweezer system according to 
any of the previous claims. 

10. A laser machining tool according to any of the previous claims. 

1 1 A method of synthesizing an output electromagnetic field u{x" f y'\ z"), comprising 
the steps of 

phase modulating an input electromagnetic field by phasor values e** ( * ,y) , 

Fourier or Fresnel transforming the phase modulated electromagnetic field, 

filtering the Fourier or Fresnel transformed electromagnetic radiation by 

in a region of spatial frequencies comprising DC in the Fourier or Fresnei 
plane 

phase shifting with a predetermined phase shift value 6 the modulated 
electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant S, and 

in a region of remaining spatial frequencies in the Fourier or Fresnei plane, 

multiplying the amplitude of the modulated electromagnetic radiation 
with a constant A 

forming an electromagnetic field o(x\ /) by Fourier or Fresnel transforming the phase 
shifted Fourier or Fresnel transformed electromagnetic field, and 

phase modulating the electromagnetic field o<x\ Y) Into the output electromagnetic 
field o(x', y')^ (K,y,) } propagating as the desired output electromagnetic field 

12. A method according to claim 1 1 , further comprising the steps of 

dividing the electromagnetic field o(x , ,y 1 ) into pixels in accordance with the disposition 
of resolution elements (x, y) of a first phase modifying element having 
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a plurality of Individual resolution elements (x, y), each resolution element 
(x, y) modulating the phase of electromagnetic radiation incident upon it with 
a predetermined phasor value e ,6w) , 

calculating the phasor values e Wxy) of the phase modifying element and the phase 
5 shift value 8 substantially in accordance with 

wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction region, 

10 B is an optional amplitude modulation of the spatial phase filter in the zero-order 
diffraction region, 

a - |a| exp (/^ ) is ths average of the phasors e* w of the resolution elements of the 
phase modifying element, and 

^ = and 
15 K m 1 - </„ (1 .22^77) * wherein 

Jo fe the zero-order Bessel function, and 

;/ relates the radius Rt of the zero-order filtering region to the radius R s of the main- 
lobe of the Airy function of the input aperture, ij^R^I^^ (0.61)" 1 ArAf t , , 

selecting, for each resolution element, one of two phasor values which represent a 
20 particular grey level, and 

supplying the selected phasor values e 1 *^ to the respective resolution elements 
(x, y) of the first phase modifying element, and 

supplying selected phasor values e W) to respective resolution elements (x\ y") of a 
second phase modifying element having a plurality of Individual resolution elements 
25 (x' ( y')> each resolution element (x\ y 3 ) modulating the phase of electromagnetic 

radiation incident upon it with the respective phasor value e^W for generation of the 
output field o(x\yy { * * 
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